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INTRODUCTION 


A.  MOTIVATION 

With  the  expanding  application  of  lasers  by  the 
military^  the  challenge  of  predicting  atmospheric  effects 
on  these  systems  has  attracted  increasing  interest.  The 
dynamics  of  an  ever-changing  atmosphere  have  profound 
effects  on  the  performance  of  high-energy  laser  (HEL) 
systems  in  the  form  of  transmission  loss^  scattering^ 
absorption^  obstruction^  and  optical  turbulence 

degradation.  Whether  applied  to  missile  defense^  laser 

communications^  or  weapon's  system  targeting^  optical 
turbulence  degradation  will  act  to  limit  the  systems 
effectiveness  and  range. 

The  capability  to  predict  HEL' s  performance  under 
specific  atmospheric  conditions  could  dramatically  enhance 
the  system  viability.  To  do  this^  higher  quality  modeling 
and  forecasts  of  the  atmosphere's  temporal  variability^ 
with  emphasis  on  optical  turbulence^  are  required.  In 

return^  military  meteorologists  will  be  given  the  ability 
to  provide  quantitative  forecasts  of  HEL' s  maximum 
effective  range  and  dwell  time.  This  thesis  will  explore 
application  of  the  U.S.  Navy's  Coupled  Ocean/Atmospheric 
Mesoscale  Prediction  Systems  (COAMPS)  model  as  a  means  to 
meet  this  requirement  and  ultimately  increase  environmental 
support  to  the  warrior. 

B .  FOCUS 

An  abundance  of  research  has  been  conducted  on  the 
application  of  high-energy  lasers  (HEL)  for  the  U.S.  armed 
forces.  In  recent  years^  significant  emphasis  has  been 
placed  on  the  negative  impacts  of  atmospheric  turbulence 
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and  its  effects  on  laser-based  weapons  and  communications 
systems.  With  key  focus  areas  being  space^  air^  and 

surface  based  elements^  each  warrants  special  attention 
when  considering  the  source  level  of  initial  transmission^ 
path  taken^  and  target  location.  Although  much  of  the 
information  provided  in  the  following  sections  is  general 

enough  to  address  the  effects  of  turbulence  in  all  of  these 
elements^  the  direction  of  this  study  approachs  the  optical 
turbulence  problem  from  the  space-based  laser  perspective 
with  emphasis  on  the  modeling  and  forecasting  of  optical 
turbulence  in  the  upper  tropospheric/lower  stratospheric 
region.  This  is  accomplished  by  comparing  a  modified 

version  of  the  COAMPS  2.0.16  mesoscale  model  output  to 
thermosonde  measurements  taken  over  Vandenberg  AFB  during 

HEL  test  periods  in  October  2001  and  March  2002. 
Specifically^  turbulence  parameterizations  are  evaluated 
and  then  adjusted  for  potential  improvements.  To  better 
understand  the  complexity  of  the  issue^  a  background  in 
atmospheric  turbulence  and  modeling  is  provided  in  the 

following  section . 
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II.  BACKGROUND 


A.  THE  CHALLENGE  OF  MODELING 

The  question  we  must  always  be  asking  is^"  What  is  the 
required  resolution  for  the  given  forecast?"  In 

operational  forecasting^  we  encounter  a  direct  conflict 
between  simulating  atmospheric  physics  accurately  and 
producing  timely  forecast  products;  this  is  addressed  by 
running  the  numerical  models  at  coarser  resolutions^  in 
both  the  vertical  and  horizontal^  as  well  as  focusing  the 
model  on  a  limited  area.  To  adequately  represent  the 

atmosphere  at  scales  smaller  than  the  resolution  being 
modeled^  sub-grid  physics  schemes^  or  parameterizations^ 
are  required.  When  running  the  model  over  a  limited  area^ 
mesoscale  models  become  extremely  sensitive  to  lateral 

boundary  conditions  in  addition  to  their  sensitivity  to 
initial  conditions^  found  in  all  models.  Generally^ 

initial  conditions  are  taken  from  the  coarser  horizontal^ 
vertical^  and  temporal  resolutions  of  synoptic  scale 
models^  thus  requiring  interpolation  into  the  mesoscale. 
With  interpolation  comes  the  potential  for  introducing 
errors  at  each  time  step  that  could  have  even  greater 
impacts  on  the  smaller  scaled  features  that  were 

parameterized.  By  understanding  the  processes^ 

sensitivities^  and  tendencies  of  a  models  these  errors  can 
be  mitigated. 

It  is  also  important  to  note  that  using  a  model  to 

research  a  phenomenon  in  retrospect  is  approached  much 
differently  than  using  a  model  to  predict  the  same 
phenomenon  in  an  operational  environment.  As  explained  by 
Warner  et  al .  (1997)^  a  subjective  estimate  of  the 


3 


importance  of  error  can  be  based  on  the  degree  to  which  the 
meteorology  in  a  given  case  is  dominated  by  initial 
conditions^  local  forcing  and  the  propagation  of  the 
features  outside  the  area  of  interest.  When  applied  to  a 
single  meteorological  case  as  is  typically  done  in 
research^  the  estimate  is  based  upon  the  prevailing 
meteorological  condition.  In  an  operational  case^ 
decisions  are  made  based  on  the  worst-case  estimates 
encompassing  an  ensemble  of  meteorological  situations. 

Addressing  our  specific  interest  in  optical 
turbulence^  determining  the  threshold  for  optimum 
timeliness  at  the  minimum  resolution  depends  on  accurate 
parameterization  of  the  turbulent  processes.  Paying 
particularly  close  attention  to  variables  within  the 
prognostic  Turbulent  Kinetic  Energy  Equation^  and  their 
inter-relationships^  will  help  us  in  this  endeavor. 


B .  TURBULENCE 

The  atmosphere  is  a  fluid  in  turbulent  motion  that 
continually  undergoes  fluctuations  in  density  from  various 
influences.  Turbulence  is^  in  the  most  general  sense^ 
random  irregularities  in  the  atmospheric  flow.  The  major 
origins  of  these  irregularities  arise  from:  (1)  mechanical 
mixing  induced  by  orographic  effects;  (2)  thermal 
interaction  associated  with  the  vertical  currents  caused  by 
surface  heating  in  the  boundary  layer;  (3)  thunderstorm 
induced  vertical  motions;  (4)  mountain  waves;  and  (5)  wind 
shear  effects.  The  latter  is  most  relevant  to  this  paper 
as  it  specifically  addresses  the  turbulence  caused  by  the 
shearing  of  horizontal  winds  in  the  vertical  direction 
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which  occur  along  air-mass  boundaries^  temperature 
inversions^  and  near  jet  streams  (Lee  and  Crane^  1994) . 

Atmospheric  turbulence  leads  to  the  development  of 
optical  turbulence  in  the  free  atmosphere  through  the 
energy  cascade  process^  optical  turbulence  occurring  on  a 
similar  scale  to  clear  air  turbulence  (CAT) .  Though  not 
identical  to  CAT^  the  temperature  and  velocity  fluctuations 
are  so  closely  linked  that  indicators  for  CAT  can  be  used 
to  identify  regions  where  optical  turbulence  would  likely 
be  found.  As  described  by  Beniston  (1998)^  turbulent 

flows  are  always  dissipative^  since  viscous  shear  stresses 
perform  deformation  work  that  increases  the  internal  energy 
of  the  fluid  at  the  expense  of  the  kinetic  energy  of  the 
turbulence.  These  characteristics  have  a  profound  effect 
on  the  dispersion  of  heat^  moisture^  and  momentum.  By 

definition^  optical  turbulence  manifests  itself  as  the 
temporal  and  spatial  fluctuations  of  the  index  of 
refraction  {n)  .  These  fluctuations  can  cause  the  HEL  beam 
to  jitter  on  the  target  as  well  as  decrease  the  energy 
density  on  target. 

A  major  source  region  of  clear  air  turbulence 
(approximately  60%) ^  and  subsequently  optical  turbulence 
through  similarity^  lies  on  the  fringes  of  the  jet  stream 
in  the  vicinity  of  upper  level  frontal  zones.  Here^  the 
wind  shear  and  thermal  gradients  act  to  produce  areas 
conducive  to  the  production  of  turbulence.  These  regions 
include  (Meteorological  Office  College^  1996) : 

1.  the  cold  side  of  the  jet^  near  and  below  the 

core; 

2.  the  warm  side  of  the  jet^  above  the  core; 
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3.  marked  curvature  and  diffluent  areas  associated 
with  jet  exit  regions; 

4.  the  confluent  and  diffluent  regions  of  two  jets; 

5.  in  the  vicinity  of  sharp  upper  troughs; 

6.  in  the  vicinity  of  the  warm  side  of  sharp  ridges; 

7.  where  one  jet  undercuts  another;  and 

8.  in  regions  of  tropopause  height  fluctuations. 

Typical  scenarios  for  development  of  turbulence  also 
include  the  amplification  of  gravity  waves  as  strong  winds 
cross  mountain  ranges  and  propagate  vertically. 
Atmospheric  waves  form  above  and  downwind  of  the  barrier  in 
a  stable  environment  as  strong  winds  travel  across 

topographic  barriers.  In  certain  thermal  conditions^  the 

resulting  waves  can  break  down  into  smaller  scale 
turbulence  (Ellrod^  2004)  .  Another  prolific  producer  of 
turbulence  is  found  in  the  region  where  a  cut-off  low 
pressure  systems  is  located  to  the  south  of  a  blocking  high 
pressure  cell  (Rosendal^  2004)  .  Additionally^  the  stably 
stratified  conditions  in  the  lower  stratosphere  can  support 
sheets  of  strong  optical  turbulence  (Mahalov  and 

Nicolaenko^  2004) . 

When  a  HEL  encounters  the  irregularities  associated 
with  turbulence^  the  optical  wave  fronts  are  subject  to 

distortion.  This  distortion  is  the  result  of  variations 
in  the  refractive  index  caused  by  the  density  fluctuation. 
These  density  fluctuations  are  produced  by  the  presence  of 
pressure  and  temperature  gradients.  The  resultant  density 


changes  occur  because  density  is  proportional  to  the  ratio 

of  pressure  and  absolute  temperature  via  the  ideal  gas  law. 
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Kyle  (1993)  defines  the  index  of  refraction  {n)  for 

wavelength  in  terms  of  pressure  (P)  in  millibars  and 

temperature  (T)  in  Kelvin^  that  are  within  or  near  the 

visible  range^  like  those  used  in  HEL  applications^  as: 

„-l  =  77.6xl0-"(l  +  3i^^)f  .  (1) 

T 

The  effects  of  turbulence  induced  pressure  variations 
are  negligible  when  compared  to  temperature  and  humidity 
fluctuations^  temperature  being  more  important  in  the  upper 
atmosphere.  While  considering  fluctuations  in  the  index  of 
refraction^  the  wavelength  can  be  ignored  because 
dependence  on  the  averaged  portion  of  the  index  becomes 
more  important  (Wolfe  and  Zissis^  1978)  .  Thus  n 

fluctuations  depend  only  on  temperature  so  that  we  can 
express  the  index  fluctuations  with  respect  to  temperature 
variations  by: 

dn  _  il9xlO-^)P 
dT  ' 

where  P  is  the  local  air  pressure  in  millibars  (mb)  and  T 
is  the  absolute  temperature  in  Kelvin. 

By  incorporating  Kolmogorov' s  theory  of  turbulence 
(Kolmogorov^  1941)^  it  is  assumed  that  the  kinetic  energy 
associated  with  larger  eddies  is  redistributed  without  loss 
to  successively  smaller  eddies  until  ultimately  being 
dissipated  by  viscosity.  Furthermore^  Kolmogorov 

determined  that  the  structure  function  for  temperature 
fluctuations  obey  the  isotropic  two-thirds  power  law^ 

Dy.(r)  =  gr^  ,  (3) 
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where  Dj.{r)  represents  the  structure  function^  is  the 
optical  turbulence  temperature  structure  parameter^  and  ris 
the  separation  distance  between  two  measurements^  in  this 
case  temperature^  that  fall  between  the  inner  and  outer 
scales  of  turbulence.  This  leads  to  the  proportionality 
relationship  between  the  optical  turbulence  structure 
parameter  as  a  function  of  the  index  of  refraction  ( )  and 

the  optical  turbulence  temperature  structure  parameter  ( ) 
given  by : 


c: 


(79xlO"*')P 


n2 


Cl 


(4) 


This  relationship  is  arrived  at  by  applying  Eg.  (3)  to  both 
cl  and  cl  where ^ 


(5) 


Cl 


(6) 


with  and  representing 

average  of  the  index  of  refraction^ 
temperature^  7^  and  ,  between  two  points, 
and  (6)  can  then  be  combined  as 


the  ensemble 

and  ,  and 
Equations  (5) 


(7) 


Despite  the  considerable  complexity  in  the  dynamics  of 
turbulence^  equations  (5)  and  (6)  give  parameters  which 
allow  us  to  identify  it  in  the  free  atmosphere.  In  situ 


measurements  of  can  be  made  with  aircraft  and  balloon- 

borne  temperature  sensors  while  Cl  measurements  can  be 
obtained  with  radar^  lidar^  and  acoustic  sounders. 

C .  NUMERICAL  MODELING 

To  adequately  model  the  dynamics  and  thermodynamics  at 
the  scales  consistent  with  turbulence^  the  fundamental 
physics  and  mathematical  laws  that  govern  the  atmosphere 
must  be  applied  in  a  manner  that  will  yield  meaningful 
results  while  maintaining  computational  efficiency.  In  a 
mesoscale  models  this  is  achieved  through  treating  the 
phenomena  as  a  parameterized  sub-gridscale  process. 
Parameterizat ions  are  used  when: 

•  the  process  being  considered  is  smaller  than  the 
characteristic  grid  of  the  models  horizontally  or 
vertically; 

•  the  physics  describing  the  phenomenon  would  act 
to  decrease  computational  efficiency  if  they  were 
computed  explicitly; 

•  or  the  scales  of  fundamental  forcing  related  to 

the  phenomenon  are  not  well  observed  or 

understood . 

Focusing  attention  on  turbulent  flow^  its 
characteristic  randomness  makes  a  deterministic  approach  to 
the  problem  implausible;  therefore^  turbulence  is  treated 
with  statistical  methods.  Fortunately^  turbulent  energy 
cascades  from  larger  eddies  to  smaller  eddies^  allowing  the 
turbulence-forecasting  to  be  an  achievable  goal  if  the 
larger-scale  numerical  prediction  is  sufficiently  accurate. 

For  this  study^  optical  turbulence  predictions  were 
performed  by  the  U.S.  Navy's  Coupled  Ocean  Atmosphere 
Mesoscale  Prediction  System  (COAMPS)  the  U.S.  Navy 
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mesoscale  model  and  compared  to  observations  at  Vandenberg 
AFB  during  October  2001  and  March  2002. 


D .  COAMPS 

A  detailed  description  of  COAMPS  is  provided  by  Hodur 
(1997)  so  that  only  a  summary  of  the  governing  equations 
and  a  brief  overview  of  some  of  its  numerical  schemes  will 
be  provided  here.  COAMPS  is  a  meso-y^  scale  model  that 
utilizes  nonhydrostatic  equations  (Klemp  and  Wilhelmson^ 
1978)  to  predict  idealized  or  real-data  simulations.  The 
importance  of  the  non-hydrostatic  condition  is  that  it 
allows  the  model  to  deal  with  convection  and  smaller-scale 
topographic  features  so  that  vertical  accelerations  are  not 
ignored.  This  is  particularly  important  for  this  study 
since  the  study  requires  modeling  at  resolutions  where  the 
hydrostatic  approximation  breaks  down  and  the  area  of 
interest  is  in  the  vicinity  of  a  coastal  mountain  range. 

COAMPS  utilize  the  Arakawa  C  grid  structure  and  has 
the  ability  to  nest  smaller  meshes  at  a  3:1  ratio.  Static 
inner  grids  can  be  arbitrarily  specified  within  the 
confines  of  the  next  coarser  grid.  in  this  case^  an  81  km 
grid  spacing  was  defined  in  the  outer  domain  with  an 
innermost  domain  grid  spacing  of  3  km.  Vertical  levels 
are  defined  in  the  sigma-z  coordinate  system  (Gal-Chen  and 
Somerville^  1975)  .  The  vertical  deminsion  is  resolved  in 
47  non-standard  levels  from  sea  level  up  to  20  km^  11  of 
which  in  the  lowest  2000  m  of  the  troposphere.  The 
vertical  resolution  in  the  upper-troposphere  and  lower- 
stratosphere  is  500  m. 
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nested  grids  over  Vandenberg  AFB . 


The  model  was  initialized  with  the  Nuss  and  Titley 
(1994)  two-dimensional  multiquadric  interpolation  scheme 
using  available  observations  blended  with  NCEP's  ETA  22  km 
horizontal  resolution  and  39  vertical  p-levels  grid  fields. 
The  verification  observations  were  not  included  in  the 
initialization  process.  All  model  simulations  begin  with 
a  '''cold"  starts  meaning  that  neither  clouds  nor  three- 
dimensional  mesoscale  structures  of  wind  or  thermal  fields^ 
were  passed  on  to  the  initialization  from  model  first  guess 
fields . 
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E.  MODELING  OPTICAL  TURBULENCE 

Since  optical  turbulence  is  not  explicitly  addressed 
within  the  COAMPS  model  physics^  it  must  be  predicted  using 
prognostic  variables  from  the  model.  As  mentioned  before^ 
the  full  set  of  prognostic  turbulence  equations  would  be 
too  costly  on  processing  time  to  be  viable  for  operational 
purposes.  To  remedy  this^  COAMPS  incorporates  Mellor  and 


Yamada  (MY)  level 

2 . 5  parameterization 

to  estimate 

the 

vertical  turbulent 

flux 

terms  u'w' ,  v' 

w' ,  and  0'w' 

as 

explained  in  Mellor 

(1973) 

,  Mellor  and  Yamada  (1974), 

and 

Mellor  and  Yamada 

(1982) . 

Even  though 

calculations 

for 

the  required  prognostic  variables  are  directly  affected  by 
closing  the  set  of  prognostic  turbulence  equations  with 
parameterization^  it  is  a  computational  necessity.  In 
this  parameterization  the  horizontal  wind  components  (U,V) 
and  potential  temperature  (0)  define  the  vertical  turbulent 
fluxes  as^ 


u  w 


-K^ 


dU 

dz' 


(8) 


V  w 


(9) 


d0 


(10) 


where  both. 


K  =  SLF 


1/2 


(11) 


and 


K,-S,LE 


1/2 


(12) 
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define  the  turbulent  momentum  and  thermal  dif fusivities^ 


respectively.  L  is  a  mixing  length  scale  and  E  is  the 
turbulent  kinetic  energy  (TKE)  such  that^ 


E  =  ^(u^  +  v"  +  w^). 


(13) 


Mellor  and  Yamada  define  as  the  momentum  stability 
function  and  Sh  as  the  heat  stability  function.  Both  are 
functions  of  gradient  Richardson  number. 


E  = 


g 

1 

du 

ydzj 


+ 


dv 

ydzj 


(14) 


The  gradient  Richardson  number  represents  the  ratio  of 
the  energy  extracted  by  buoyancy  forces  to  the  energy 
gained  from  the  shear  forces^  giving  quantitative 
information  on  the  relation  between  the  stabilizing  effect 
of  buoyancy  and  the  destabilizing  effect  of  velocity  shear. 
Once  the  ratio  exceeds  a  critical  threshold  of  0.25^ 
turbulence  can  no  longer  be  generated.  This  particular 
study  will  focus  on  examining  stable  atmospheric 
conditions . 


Given  the  proportional  relationship  between  the 
optical  turbulence  parameter  ( )  and  thermal  structure 

parameter  ( )  seen  in  (3)^  the  problem  can  be  approached 
by  solving  for  .  In  a  stable  atmosphere  where 

depends  on  the  atmospheric  stratification^  the  thermal 
structure  parameter  can  be  expressed  in  terms  of  the  bulk 
potential  temperature  {6)  as  explained  by  Bougeault  et  al . 
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(1995) . 


The  resultant 


in  terms  of  the  mechanical 


mixing  length  scale  {L^)  and  the  dissipation  length  scale 
(4)  is 


4  - 1  4Z  { — 

V  OZ 


This  was  further  simplified  in  Walters  and  Miller  (1999)  in 
terms  of  a  special  mixing  length  scale  (4^  that 


To  date^  the  TKE  and  TKE-Free  Methods  show  great 
potential  for  aiding  us  in  the  forecasting  of  optical 
turbulence.  The  TKE  Method^  as  described  in  Walters  and 
Miller  (1999)^  will  be  examined  in  this  thesis.  This 

method  uses  three-dimensional  predictions  of  potential 
temperature  and  turbulent  kinetic  energy^  one  of  the 
prognostic  variables  affected  by  the  Mellor  and  Yamada  (MY) 
parameterization  scheme.  The  TKE-Free  Method^  which  is 

important  to  mention  for  future  investigation^  uses 
predictions  of  potential  temperature  and  winds.  Each 

method  calculates  4  ^  different  manner^  as  explained 

below . 

1 .  TKE  Method 

To  achieve  realistic  turbulent  kinetic  energies  and 

temperature  fluctuations  in  the  stable  atmosphere  an 

appropriate  thermal  to  eddy  diffusivity  (Kh/K^)  dependence 

must  be  incorporated  into  the  turbulent  kinetic  energy 

tendency  equation.  Although  the  default  MY  eddy 
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meet 


this 


diffusivity  parameterization  attempts  to 
requirement,  known  deficiencies  exist  and  are  addressed 
later  in  this  section.  The  COAMPS  prognostic  equation  for 
TKE  (Hodur  1997)  is 


dt 


dE  dE  dE  ,, 

- + 

OX  oy  oz 


\dzj 


ydzj 


g  ^0 

e  dz 


-X^—  +  D+  KX^E  -  £. 


(17) 


(a) 


(b)  (c)  (d)  (e)  (f) 


Components  of  (17)  are: 

(a)  three-dimensional  advection; 

(b)  shear  production; 

(c)  buoyancy  production; 

(d)  sub-grid  scale  mixing; 

(d)  numerical  diffusion; 

(e)  and  dissipation. 

The  thermal  mixing  length  applied  in  the  TKE  Method 
uses  the  Deardorff  (1980)  definition  where 


L  =  0.76 


(18) 


With  regards  to  eddy  diffusivity  parameterization^ 
Walters  (1995)  and  Walters  and  Bradford  (1997)  noted  that 
when  comparing  observed  optical  turbulence  TKE  profiles 
with  modeled  TKE  values^  more  realistic  results  where 
achieved  by  using  the  Kondo  et  al .  (1978)  polynomial^ 
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_  «1 

i  ■ 

a,R.+ - 

'  \  +  a^R^ 

Miller  et  al ,  (2001)  showed  that  the  simulated 
turbulent  processes  are  quite  sensitive  to  the  K^jK^vs.  R. 
relationship.  Generally^  too  aggressive  of  a 

parameterization  can  exaggerate  the  turbulent  processes 
while  a  too  conservative  approach  can  dampen  it  out 
completely . 


2.  TKE-Free  Method 

This  method  represents  a  lower  order  closure  of  the 
vertical  turbulent  fluxes  compared  to  the  TKE  Method  and^ 
as  stated  before^  uses  predictions  of  potential  temperature 
and  winds.  It  estimates  the  mixing  length  scale  required 
in  Eq.  (16)  based  on  observations  presented  by  Tjernstrom 
(1993)  so  that 


(l  +  l5Ri^l  +  5R. 


where  Ri  is  defined  in  Eq.  (14) . 
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III.  ORGANIZATIONAL  APPROACH 


Initially^  a  synoptic  re-analysis  of  the  Eastern 
Pacific  is  provided  for  both  the  October  2001  and  March 
2002.  Emphasis  will  be  placed  on  those  synoptic  scale 

features  that  would  most  likely  produce  optical  turbulence 
over  the  area  of  interest^  Vandenberg  AFB . 

To  gauge  of  the  model’s  initial  performance^  COAMPS' s 
horizontal  model  output  will  be  directly  compared  to  the 
synoptic  analysis  from  National  Centers  for  Environmental 
Prediction/Marine  Prediction  Center  (NCEP/MPC)  at  OOZ  and 
12Z^  or  closest  available  analysis  if  those  times  were  not 
available.  This  will  validate  the  model's  accuracy  and 

expose  inconsistencies  that  may  have  been  inherited  by 
errors  in  lateral  boundary  conditions  or  by  poor 
initialization.  Additionally^  rawinsonde  data  will  be 

used  to  compare  vertical  depictions  of  the  model  output  at 
the  closest  relevant  grid-point.  These  steps  will  reveal 
inconsistencies  that  may  be  present  in  the  column  above  the 
point  of  interest.  Next^  the  TKE  Method  will  be  applied 
to  3  km  CCAMPS  runs  of  the  Cctober  2001  case  to  produce 
synthetic  vertical  profiles.  This  will  be  accomplished  by 
creating  a  synthetic  balloon  launch  within  the  model  to 
produce  a  trajectory  comparable  to  that  of  an  actual 
profile.  These  profiles  will  in  turn  be  compared  to 
corresponding  ''real  world"  balloon  launched  thermosonde 
data  and  qualitatively  analyzed  in  detail.  Additionally^ 
these  results  will  be  quantitatively  compared  by 
determining  the  isoplanatic  angle^  defined  in  CH .  V  Sec. 
C.^  for  each  sonde  to  produce  a  scatter  plot.  Results 
will  subsequently  be  used  to  make  adjustments  to  measure 
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the  improvement  of  new  parameterizations  in  outlying  model 
runs  selected  from  the  October  2001  and  March  2002  cases. 
The  "Default"  and  "Adjusted"  COAMPS  model  runs  for  October 
2001  and  March  2002  will  then  be  quantitatively  compared  to 
thermosonde  data  in  the  same  manner  as  described  for  the 
October  2001  case. 
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IV.  SYNOPTIC  ANALYSIS 


A.  OCTOBER  2001 

For  this  period  of  study^  the  synoptic  conditions  for 
October  18^^  through  25^^  are  derived  from  the  National 
Climatic  Data  Center's  (NCDC)  surface  analysis^ 

supplemented  by  Unysis  composite  GOES  infrared  satellite 
imagery  and  300  mb  AVN  jet  level  wind  analysis  from  San 
Francisco  State  University  Meteorology  Department. 
Conditions  conducive  to  the  production  of  optical 
turbulence  will  be  highlighted  to  help  focus  mesoscale 
analysis  on  those  periods  likely  to  have  affected  the 
experiment  over  Vandenberg  AFB .  Relevant  charts  and 

imagery  will  be  shown  in  the  following  synoptic  summary 
while  the  complete  compliment  of  charts  are  located  in 
Appendix  A  with  satellite  imagery  in  Appendix  B. 

Beginning  on  October  18th^  the  synoptic  picture  is 
dominated  surface  high-pressure  system  about  800  nautical 
miles  (nm)  west  of  Oregon  as  can  be  seen  in  Figure  2 . 
This  persistent  high  remains  a  synoptic  forcing  feature 
throughout  the  period.  A  semi-permanent  thermal  low  in 
southern  California  is  present  in  the  beginning  of  the 
period  with  an  associated  inverted  trough  extending  up 
though  the  California  Central  Valley^  affecting  weather 
over  the  area  of  interest.  The  300  mb  jet  level  winds^  as 
depicted  in  the  AVN  analysis  in  Figure  3^  show  the  Polar 
Front  Jet  (PFJ)  north  of  the  surface  high  arcing  from  the 
Gulf  of  Alaska  into  British  Columbia. 
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By  12Z  on  the  18th^  a  cyclonic  circulation  off  the 
coast  of  California^  indicative  of  an  upper  level  cut-off 
low^  becomes  apparent  in  the  jet  level  winds  depicted  in 
the  300  mb  AVN  analysis.  Satellite  imagery  from  the  GOES 
composite  indicates  that  a  well-developed  circulation  had 
developed  from  12Z  on  the  18th  to  OOZ  on 
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the  19th  with 


high-level  clouds  exhibiting  the  characteristic  ’comma’ 
signature  associated  with  cyclonic  circulation  (Figure  4). 
On  the  surface  analysis^  this  upper  level  feature  is 
manifested  as  a  stationary  trough  (Figure  5) . 


Figure  4.  GOES  OOZ  Composite  IR,  19OCT01. 


Figure  5.  NCEP/MPC  12Z  Surface  Analysis,  19OCT01. 
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As  the  Northeast  Pacific  High  begins  to  weaken 
slightly  on  the  19th^  the  300  mb  winds  on  the  southeast 
quadrant  of  the  cut-off  low  begin  to  intensify.  This 
intensification  is  the  result  of  two  different  processes 
contributing  to  its  development.  Firsts  cold  air  advected 
by  the  Northeast  Pacific  high  is  being  wrapped  up  below  the 
upper-level  cutoff  low^  helping  to  deepen  it.  Second^ 
interaction  with  the  Subtropical  Jet  (STJ)  produces  speed 
convergence  along  the  southern  edge  of  the  cut-off  low. 
These  two  interactions  cause  a  jet  maximum  to  form  near  30 
N  135  W.  By  the  21^^^  the  jet  extends  to  the  central  coast 
of  California  and  begins  to  propagate  eastward  as  the 
Northeast  Pacific  high  starts  rebuilding.  The  jet  edge  is 
clearly  visible  crossing  central  California  on  the  Figure  6 
satellite  image. 


Figure  6.  GOES  12Z  Composite  IR,  21OCT01. 

Continuing  through  the  period^  the  jet  max  moves 
eastward  until  it  merges  with  the  PFJ  on  the  23rd.  At  the 
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end  of  the  period^  the  northeast  Pacific  high  moves 
southwest  of  its  original  position  from  the  18th  and 
reestablishes  itself  as  the  dominant  weather  influence  over 
the  California  coast.  The  PFJ  has  similarly  reoriented 
itself  northwesterly  across  Oregon  and  Washington  state. 

This  STJ  core  is  of  particular  interest  because  it  is 
expected  to  be  a  prolific  producer  of  turbulence  in  the 
upper  levels  over  our  area  of  interest.  Similarly^  since 
the  beginning  of  the  period  was  not  affected  by  a  strong 
turbulence  producing  feature^  this  event  should  provide  a 
an  opportunity  to  test  the  forecast  capability^  both 
strengths  and  weaknesses^  of  mesoscale  COAMPS  analysis. 

B.  MARCH  2002 

This  second  period  covers  synoptic  features  over  the 
Eastern  Pacific  from  March  1^^  through  7^^^  2002.  Surface 
analysis  charts  used  in  this  review  are  the  eastern  segment 
of  NCEP's  Pacific  surface  analysis.  Both  NCEP  analysis 
and  infrared  4  km  GOES  satellite  half  disk  imagery  were 
obtained  from  archives  of  the  Naval  Postgraduate  School 
(NPS)  Meteorology  Department.  Jet  level  winds  were  taken 
from  the  300  mb  AVN  analysis  archived  by  San  Francisco 
State  University  Meteorology  Department. 

On  March  1^  the  surface  analysis  shows  a  high-pressure 
system  just  west  of  the  Washington  coast  and  a  stationary 
trough  between  Hawaii  and  California.  Over  southern 
California^  a  stationary  low  is  anchoring  an  inverted 
trough  that  extends  to  the  north  (Figure  7)  .  This  surface 
pattern  persists  until  March  4^^. 
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Figure  7.  NCEP/MPC  12Z  Surface  Analysis,  01MAR02. 


Figure  8.  AVN  12Z  300mb  Jet  Stream  Analysis,  01MAR02. 

Aloft,  Figure  8  shows  that  the  PFJ  is  highly 
meridional,  becoming  northerly  over  Washington  and  Oregon 
before  deviating  sharply  to  the  East  over  northern 
California.  The  STJ  is  zonal  along  the  Baja  Peninsula  and 
northern  Mexico.  Continuing  through  the  period,  the  PFJ 
extends  further  south  along  the  California  coast, 
eventually  merging  with  the  STJ  on  the  2^^,  then  becoming 


southwesterly  over  the  middle  of  the  United  States  (Figure 
9)  .  IR  satellite  imagery  depicts  this  feature  with  clear 
skies  to  the  north  of  the  jet  and  scattered  cloudiness  on 
the  warm  side  (Figure  10),  as  seen  in  imagery  from  the 
to  the  5’^'^. 


Figure  9.  AVN  12Z  SOOmb  Jet  Stream  Analysis,  02MAR02 . 


Figure  10.  GOES  WEST  18Z  IR,  02MAR02 . 
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On  the  4^^  of  March^  the  pattern  begins  to  change  as  a 
short  wave  interacts  with  the  stationary  trough  that  was 
located  near  140W^  resulting  in  development  of  the  low  and 
f rontogenesis .  This  new  baroclinic  system  with  its 

associated  fronts  begins  tracking  northeast.  Due  to 

interaction  with  the  STJ^  a  jet  maxima  forms  in  the 
southeastern  quadrant  of  the  cyclonic  circulation  as  it 
intensifies^  extending  the  jet  northeastward  as  with  the 
frontal  system  propagates  (Figure  11).  This  elongation  of 
the  jet  continues  until  it  crosses  the  California  coast  on 
5^^  and  persists  until  mid-period  (12Z)  on  the  7^^. 


Figure  11.  AVN  18Z  300mb  Jet  Stream  Analysis,  06MAR02. 

The  STJ  interacts  with  the  PFJ  on  the  8^^  and  the  two 
jets  again  combine  to  form  a  300  mb  trough  that  quickly 
propagates  eastward,  setting  the  stage  for  another  low- 
pressure  system  to  move  into  the  California  coast  mid¬ 
period  on  the  9^^. 
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In  summary^  the  PFJ  placement  on  the  1^^  and  2^^  of 
March  are  expected  to  produce  optical  turbulence  aloft. 
From  the  3^^  to  late  on  the  5^^^  conditions  aloft  are  less 
conducive  to  produce  extreme  gradients  in  the  index  of 
refraction.  During  the  second  half  of  the  study^  when  the 
STJ  crosses  the  West  Coasts  we  would  once  again  expect 
higher  values  of  optical  turbulence^  including  the 
interaction  that  takes  place  between  the  PFJ  and  STJ  on  the 
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V. 


MODEL  ANALYSIS 


A.  COAMPS  ANALYSIS 

1 .  Surface  Analysis 

The  first  analysis  that  was  performed  compared  the 
NCEP  Pacific  OOZ  and  12Z  Analysis  with  the  corresponding 
COAMPS  81  km  resolution  model  runs  to  ensure  that  the  model 
was  depicting  the  synoptic  situation  with  reasonable 
accuracy.  Furthermore^  both  the  previous  12  hour  forecast 
and  the  OOZ  model  analysis  valid  for  the  particular  NCEP 
surface  analysis  were  compared  and  contrasted  against  one 
another.  This  was  done  to  investigate  if  poor  turbulence 
forecast  could  be  traced  to  poor  initialization  of  the 
synoptic  conditions . 

As  can  be  seen  when  comparing  Figure  12  with  Figures 
13  and  14^  the  synoptic  NCEP  analysis  and  model  runs  match 
relatively  well  considering  the  low  resolution  depicted. 
This  example  also  shows  artifacts  that  are  common 
considering  the  forecast  time.  Specifically^  it  can  be 
seen  in  Figure  13  that  the  12-hour  forecast  shows  increased 
detail  in  isobar  structure  but  tends  to  make  the  high- 
pressure  center  slightly  higher  than  the  NCEP  analysis. 
Also  commonly  encountered  as  the  forecast  progresses^ 
distortion  of  isobaric  gradients  becomes  more  clearly 
visible  near  the  lateral  boundaries  of  the  model.  The  OOZ 
model  analysis  typically  depicts  the  high-pressure  system 
in  the  Northeast  Pacific  similar  to  the  NCEP  analysis  with 
isobars  that  more  relaxed  when  directly  compared  to  the 
prior  run's  12-hour  forecast.  Additionally^  lateral 
boundaries  for  this  run  haven't  distorted  the  isobaric 
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pattern  as  seen  in  the  12-hour  forecast.  Also  significant 
is  that  the  surface  trough  at  140W  is  wising  in  the  12  hour 
forecast . 


Figure  12 . 


NCEP/MPC  Surface  Analysis^ 


12Z  19OCT01. 
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19OCT01 . 
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A  similar  review  of  COAMPS  analysis  and  12  Hour 
forecasts  showed  satisfactory  agreement  depicting  the 
synoptic  situation  from  the  19^^  to  the  25^^  of  October. 
The  sea  level  pressures  were  within  5  millibars  at  the 
high/low  pressure  centers.  The  complete  set  of  81  km 

COAMPS  runs  for  this  time  period  are  provided  in  Appendix 
C. 


2 .  Jet  Analysis 

Jet  stream  analysis  was  performed  by  comparing  the 
AVN  300  mb  model  analysis  with  COAMPS' s  previous  12-hour 
forecast  and  the  OOZ  analysis  using  isotachs  and  winds  at 
the  300  mb  level.  This  was  conducted  in  the  same  manner 
as  the  surface  analysis  just  discussed.  Figures  15^  16 

and  17  show  the  typical  scheme  used  for  this  comparison. 


Figure  15.  AVN  Jet  Stream  Analysis  12Z,19OCT01. 
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Attention  was  directed  to  the  position  of  the  jet 
stream^  orientation^  and  the  maximum  speeds  occurring 
within  the  lateral  boundaries  of  the  model  in  relation  to 
the  maximum  speeds  and  wind  directions  taken  from  the  AVN 
300  mb  analysis. 
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Typically^  the  position  of  the  jet  maxima  in  the  12- 
hour  forecasts  were  slightly  behind^  or  more  westerly  in 
position^  to  that  of  the  AVN  analysis  while  the  analysis 
typically  had  better  position^  as  should  be  expected. 
Aside  from  the  slight  differences  in  position^  the  COAMPS 
runs  were  depicting  the  300  mb  jets  reasonably.  Applicable 
AVN  analysis  are  located  in  Appendix  D  while  the  full  set 
of  COAMPS  300  mb  forecast  and  analysis  are  located  in 
Appendix  E . 


Since  the  focus  of  this  study  delves  into  the  ability 
to  forecast  turbulence  in  the  mesoscale^  particular 
attention  was  given  to  the  times  surrounding  the  21OCT01 
12Z  analysis  over  Vandenberg  AFB^  shown  in  Figure  18. 
On  this  date^  an  extension  of  the  Sub-Tropical  Jet^  a 
prolific  producer  of  turbulence^  was  affecting  the  nested 
area  of  interest. 


Figure  18.  AVN  12Z  21OCT01  Jet  Stream  Analysis. 
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3 .  COAMPS  Skew-T  Comparison  with  Rawinsondes 

Satisfied  that  the  COAMPS  model  was  correctly 
depicting  the  synoptic  scale  thus  far^  the  next  logical 
step  was  to  verify  the  mesoscale.  This  was  achieved  by 
directly  comparing  balloon-launched  rawinsonde  measurements 
in  the  skew-T  format  with  synthetically  produced  profiles 
for  the  same  area^  Vandenberg  AFB^  using  COAMPS  forecast 
from  the  3  km  domain.  Although  not  as  detailed  as  a 

thermosonde/model  comparison  to  come^  this  step  provide  a 
more  general  indication  of  whether  the  model  is  providing 
reasonable  output . 

Table  1  shows  the  balloon  launched  that  was  compared 
to  a  particular  model  run.  The  OOZ  19OCT01  balloon  launch 
was  unavailable  so^  the  earlier  launch  (18Z  18OCT01)  was 

used.  Equally  important  to  mention^  COAMPS  runs  for 

21OCT01  12Z  and  22OCT01  OOZ  valid  times  where  not  initiated 
due  to  thermosonde  failure. 


Balloon 

ttZ  ddmmmyy 

Model 

dd  ttZvtt 

18Z  18OCT01 

18  12Zvl2 

OOZ  20OCT01 

19  12Zvl2 

00Z21OCT01 

20  12Zvl2 

OOZ  23OCT01 

22  12Zvl2 

OOZ  24OCT01 

23  12Zvl2 

OOZ  25OCT01 

24  12Zvl2 

Table  1.  Balloon  launch  times  and  the  COAMPS  forecast. 

In  all  cases^  the  12-hour  forecast  from  the  model  was 
used  for  the  comparison.  The  reasons  for  using  the  12- 
hour  forecast  instead  of  the  OOZ  model  analysis  was;  1) 
model  instability  for  the  first  six  hours  after 
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initialization^  and  2)  a  consistent  surface  temperature 
discrepancy  apparent  in  the  OOZ  analysis.  Figure  19  shows 
a  typical  comparison  of  forecast  with  the  analysis  from  the 
same  valid  times  to  demonstrate  this  point.  As  can  be 
seen^  the  surface  temperature  in  the  analysis  is  10  degrees 
Celsius  greater  than^  more  in  some  cases^  the  surface 
temperature  reflected  in  the  balloon  launched  data. 
Although  an  interesting  anomaly^  because  this  takes  place 
during  a  time  of  expected  instability^  this  observation 
does  not  interfere  with  the  focus  of  this  study^ 
forecasting  turbulence  above  the  boundary  layer. 


rOK  Oct  2)  2C01  OOOOVOI2  CQMVS4  Btow-T  ivbg' 


rcK  Oct  22  2001  oooovooo  ooftmi  (vbg;' 


Figure  19.  22OCT01  12Zvl2  COAMPS  f orecast ( lef t )  and  the 

23OCT01  OOZvOO  COAMPS  analysis  (right) . 

Figures  20  and  21  show  a  typical  comparison  between 
balloon-derived  skew-T  with  the  corresponding  forecast. 
Overall^  balloon-launched  sonde  comparisons  to  COAMPS' s 
vertical  model  where  reasonable  with  the  following 
discrepancies  noted : 

1.  The  200CT01  12Z  model  run  (Figure  152)^  12  hour 

forecast^  showed  a  marked  departure  from  the 
corresponding  balloon-sonde  (OOZ  21OCT01)  at  the  500 
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mb  levels  Figure  136.  Although  the  actual 

conditions  seem  to  indicate  a  drying  at  this  levels 
the  model's  dry  layer  is  clearly  excessive. 

2.  The  22OCT01  12z  model  run's  12  hour  forecast  (Figure 
155)^  indicated  a  dryer  column  with  consistently 
lower  dew  point  temperatures  than  that  indicated  by 
the  OOZ  23OCT01  balloon^  Figure  140. 

3.  The  23OCT01  12Z  model  run's  12  forecast  (Figure  158) 
shows  a  stronger  inversion  with  less  moisture  from 
800  mb  to  300  mb  when  compared  to  the  OOZ  24OCT01 
balloon  profile^  Figure  142.  Additionally^  stronger 
and  more  northerly  winds  are  indicated  around  the 
250  mb  level  in  the  model. 

Archived  Sounding  data  from  the  University  of  Wyoming 
is  located  in  Appendix  F  with  corresponding  COAMPS  vertical 
profiles  taken  from  the  3  km  nest  located  in  Appendix  G. 
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72393  VBG  Vandenberg  Afb 
100 


200 

300 


SLAT  34.75 
SLON  -120. 
SELV  121.0 
SHOW  6.06 
LIFT  10.33 
LFTV  10.37 
SWET  76.50 
KINX  16.50 
CTOT  11.30 
VTOT  26.30 
TOTL  30.60 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  -0909 
EQTV  -9999 
LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  261.9 
LCLP  902.3 
MLTH  290.4 
MLMR  7.97 
THCK  5715. 
PWAT  23.97 


OOZ  23  Oct  2001 


University  of  Wyoming 


Figure  20.  Skew-T  derived  from  balloon  launched  data^ 


OOZ  23OCT01. 


23OCT01 . 
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4 .  COAMPS  Skew-T  Comparison  with  Thermosonde  Data 

Similar  to  the  analysis  performed  in  Section  3^  this 
step  served  as  verification  that  the  model  was  depicting 
the  environment  well  in  a  relative  sense.  A  thermosonde  is 
a  balloon-launched  sensor  that  measures  temperature 
differences  of  less  than  one  percent  of  a  degree  centigrade 
between  two  locations  in  the  atmosphere.  Thus^  the 
thermosonde  measurements  represent  extremely  high 
resolutions  (up  to  350  data  points  per  kilometer)  in  the 
vertical  while  the  model  depiction  is  coarse  (one  point  per 
500m  in  the  upper  troposphere^  47  points  throughout  the 
column)  .  It  is  also  reasonable  to  expect  that  sondes 
launched  on  the  same  day  should  show  similar  structure 
since  they  are  launched  only  hours  apart. 

To  reiterate  an  earlier  pointy  finding  the  balance 
between  resolution^  computational  cost^  and  timeliness  of 
the  product  is  key  to  operational  viability.  Since  this 
step  simply  served  as  verification^  details  will  be 
explored  in  the  next  section.  Table  2  shows  a  list  of  the 
thermosonde^  launch  times^  and  the  applicable  COAMPS  run 
that  were  compared.  The  full  set  of  thermosonde  profiles 
are  included  in  Appendix  H. 
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The  rmo  s  onde 

Launch  time  (UT) 

Model  time  (UT) 

Thermosonde 

2 

19OCT01 

0132 

18OCT01 

12Zvl4 

Thermosonde 

3 

19OCT01 

0322 

18OCT01 

12Zvl 6 

Thermosonde 

5 

200CT01 

0115 

19OCT01 

12Zvl4 

Thermosonde 

6 

200CT01 

0300 

19OCT01 

12Zvl 6 

Thermosonde 

7 

200CT01 

0444 

19OCT01 

12Zvl7 

Thermosonde 

10 

21OCT01 

0235 

200CT01 

12Zvl4 

Thermosonde 

11 

21OCT01 

0413 

200CT01 

12Zvl 6 

Thermosonde 

12 

23OCT01 

0115 

22OCT01 

12Zvl4 

Thermosonde 

13 

23OCT01 

0254 

22OCT01 

12Zvl6 

Thermosonde 

15 

24OCT01 

0114 

23OCT01 

12Zvl4 

Thermosonde 

16 

24OCT01 

0248 

23OCT01 

12Zvl4 

Thermosonde 

17 

24OCT01 

0430 

23OCT01 

12Zvl6 

Thermosonde 

18 

25OCT01 

0115 

24OCT01 

12Zvl4 

Thermosonde 

19 

25OCT01 

0253 

24OCT01 

12Zvl 6 

Thermosonde 

21 

26OCT01 

0249 

25OCT01 

12Zvl4 

Thermosonde 

22 

26OCT01 

0423 

25OCT01 

12Zvl6 

Table  2.  Thermosonde  and  corresponding  model  run  times. 
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B.  QUALITATIVE  EVALUATION  OF  C]  PROFILES  DERIVED  FROM 
THERMOSONDE  AND  COAMPS  DATA 

This  analysis  was  performed  by  overlaying  results 
from  both  the  thermosonde  data  and  the  corresponding  COAMPS 
model  forecast.  Since  the  thermosonde  data  was  at  a  much 

higher  resolution^  as  demonstrated  in  Figure  22^  the  C'l 
values  were  averaged  over  50  data  points^  roughly  400  m^  to 
reduce  noise  while  maintaining  the  optical  turbulence  trend 
for  comparison^  as  shown  in  Figure  23. 
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Figure  23.  Thermosonde  data  averaged  over  50  points. 

Because  this  paper  is  focusing  on  the  upper- 
troposphere  interactions  affecting  optical  turbulence^  this 
analysis  concentrates  on  heights  above  the  boundary  layer 
from  3000  m  to  17500  m.  The  model's  sponge  layer  extends 
from  17.5  km  to  20  km^  making  any  quantitative  or 
qualitative  comparison  in  that  region  minimally  valuable. 
The  model  balloon  ascent  assumes  ^'steady  state"  by  using  3- 
D  model  flow  field  at  a  single  model  output  time  during  its 
ascent.  Although  steady  state  is  a  reasonable 
approximation  to  the  synoptic  conditions  shown  by 
Thermosondes  21  -  22  it  may  not  as  reasonable  approximation 
to  the  dynamic  synoptic  conditions  depicted  by  Thermosondes 
2-19.  As  a  result^  this  may  be  the  source  of  some  error 
in  the  profiles. 
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Results  from  the  surface  to  3000  as  represented  by 
Thermosonde  2  in  Figure  24^  are  not  commented  on  here  since 
this  paper  is  not  focused  in  the  boundary  layer. 


Thernnosonde:vanfa002.txt  Model:c2001 10181 2F01 6.txt 


Figure  24.  Representative  comparison  of  boundary  layer. 

Trends  from  this  analysis  will  identify  model 
discrepancies^  allowing  detailed  investigation  into  the 
possible  causes.  Profiles  that  exhibit  similar  patterns 
and  maintain  values  within  an  order  of  magnitude  are 
considered  acceptable.  Profiles  that  diverge  with  values 
exceeding  an  order  of  magnitude  are  considered 
unacceptable^  thus  model  failure.  This  in  turn  will  give 
some  insight  into  possible  corrections  that  should  be 
applied . 
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1. 


Thermosonde  2 


The  most  notable  artifact  in  Figure  25  is  the  under¬ 
forecasting  of  cl  present  just  above  10  km.  The  error  of 
the  model  is  one  order  of  magnitude.  However^  when  the  peak 
in  the  thermosonde  profile  is  compared  to  the  9  km  model 
value^  forecast  error  is  less.  Further  comparison  between 
thermosonde  and  model  derived  skew-T's^  Figures  146  and 
167^  indicate  that  the  stronger  winds  were  forecast  at  a 
higher  level  than  those  measured  by  the  thermosonde. 
Possibly  relevant  as  well^  the  model's  dewpoint  temperature 
was  higher  than  the  actual  environment^  indicating  cloud 
ice  at  300  mb  or  9km  in  the  model. 
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2. 


Thermosonde  3 


Figure  26.  Cl  profile  comparison  for  19OCT01  04Z. 


In  Figure  26^  the  model  appears  to  follow  the 
thermosonde  relatively  well  with  the  exception  of  under¬ 
forecasting  turbulence  at  11  and  13  km.  The  trend  of  the 
model  does  show  an  increase  at  12  km  but  appears  to  dampen 
it  out  too  quickly  as  the  altitude  increases.  The  only 
apparent  deviation  between  corresponding  skew-T's  is  the 
higher  dewpoint  temperature  in  the  models  Figures  147  and 
168.  Again^  the  model  does  indicate  the  presence  of  cloud 
ice  just  above  9  km. 
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3. 


Thermosonde  5 


Figure  27.  Cl  profile  comparison  for  200CT01  02Z. 

The  first  deviation  in  this  case  occurs  at  an  altitude 
around  4  km^  as  seen  in  Figure  27.  Comparison  between  the 
thermosonde  and  model  vertical  profiles  point  to  a 
significant  drying  at  the  600  mb  level  that  is  not 
represented  by  the  model's  profile^  Figures  149  and  169. 
At  8  km^  the  thermosonde  records  much  lower  dewpoint 
temperature  than  what  is  modeled  while  significant  cloud 
ice  is  shown  in  the  model. 
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4. 


Thermosonde  6 


Optical  Turbulence  (Cn^) 


Figure  28.  Cl  profile  comparison  for  200CT01  04Z. 

Unlike  the  previous  cases^  Figure  28  shows  that  the 
model  has  over-forecast  the  turbulence  at  10  and  15  km 
while  under-forecasting  it  at  4  and  6  km.  At  the  under¬ 
forecasted  levels^  the  thermosonde  again  indicates  a  much 
dryer  layer  that  what  the  model  forecast^  Figures  150  and 
170.  As  for  the  over-forecasts^  higher  dewpoint 

temperatures  are  shown  in  the  model  profile  with 
indications  of  significant  cloud  ice  around  9  km^  also 
model  winds  are  slightly  stronger  from  250  mb  to  150  mb. 
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5. 


Thermosonde  7 


Figure  29.  Cl  profile  comparison  for  200CT01  05Z. 

The  analysis  for  Figure  29  is  very  similar  to  that  of 
Figure  28^  obviously  due  to  time  similarities.  The  under¬ 
forecasted  turbulence  at  3.5  km  corresponds  to  a  much  dryer 
layer  than  what  the  model  has  depicted.  Slightly 
stronger  model  winds  and  a  warmer  dewpoint  temperature 
differentiate  the  model  profile  from  the  thermosonde 
profile  relative  to  the  over-forecast  seen  at  9  km^  Figures 
151  And  171.  This  is  complemented  by  a  strong  indication 
of  cloud  ice  around  8.5  km.  No  strong  indications  of 
vertical  wind  shear  seemed  to  be  present. 


48 


6. 


Thermosonde  10 


Figure  30.  Cl  profile  comparison  for  21OCT01  02Z. 

This  profile  comparison  shown  in  Figure  30  indicates  a 
severe  under-forecast  at  13  km.  This  anomaly  corresponds 
to  an  increase  in  temperature  on  the  sounding  that 
indicates  it  is  at  the  tropopause.  Greater  vertical  wind 
shear  is  recorded  at  this  level  than  was  modeled^  Figures 
153  and  172.  The  model  once  again  depicts  a  larger 

dewpoint  temperature  than  the  thermosonde  with  cloud  ice 
present  at  8.5  km.  The  under-forecast  at  9  km  seems  to 
correspond  with  a  less  stable  gradient  in  the  model  than 
measured . 
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7. 


Thermosonde  11 


Figure  31.  Cl  profile  comparison  for  21OCT01  04Z. 

The  most  significant  deviation  of  the  model  in  Figure 
31  appears  to  be  the  over-forecast  at  10  km.  Two  peaks  in 
cloud  ice  are  modeled  at  8.5  km  and  10  km^  the  first  of 
which  being  more  pronounced  with  cloud  indicated  on 
satellite  imagery.  The  model  indicates  a  dryer  and  less 
stable  layer  than  what  was  actually  measured^  Figures  154 
and  173.  Otherwise  the  trend  seems  reasonably  depicted. 
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8. 


Thermosonde  12 


Optical  Turbulence  (Cn^) 

Figure  32.  profile  comparison  for  23OCT01  02Z. 

With  the  exception  of  the  under-forecast  that  occurs 
at  15  km^  the  model  profile  shown  in  Figure  32  appears  to 
closely  parallel  the  thermosonde  trends  but  offset  by  about 
1  km  below  that  of  the  thermosonde.  The  only  definitive 
difference  that  can  be  noted  between  model  and  measurement 
is  that  the  temperature  gradient  of  the  model  at  15  km  is 
less  stable  in  the  models  Figures  156  And  174.  No 

conclusive  information  can  be  obtained  from  the  dewpoint 
temperature  due  to  apparent  sensor  failure. 
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9. 


Thermosonde  13 


Optical  Turbulence  (Cn^) 


Figure  33.  Cl  profile  comparison  for  23OCT01  04Z. 

Two  instances  of  under-forecasting  occur  in  the  and 
above  the  stratospheric  boundary  layer  as  shown  in  Figure 
33.  The  model  depicts  both  a  less  stable  temperature 
gradient  and  higher  dewpoint  temperature  than  the 
thermosonde^  Figures  157  And  175. 
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10.  Thermosonde  15 


Figure  34.  Cl  profile  comparison  for  24OCT01  02Z. 

The  over-forecast  that  occurs  from  11  km  to  14  km  in 
Figure  34  appears  to  correlate  with  the  more  erratic 
dewpoint  temperature  structure  of  the  thermosonde  that  was 
not  represented  in  the  moister  and  more  stable  structure  of 
the  models  Figures  159  and  176. 
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11.  Thermosonde  16 


Figure  35.  Cl  profile  comparison  for  24OCT01  04Z. 

The  only  marked  deviation  of  the  model  in  the  case 
shown  in  figure  35  appears  at  8  km.  With  the  exception  of 
weaker  winds  in  the  model  at  that  levels  no  definitive 
differences  can  be  drawn  from  the  skew-T  profiles  for  this 
occurrence^  Figures  160  And  177.  No  cloud  ice  was  modeled 
in  the  profile. 
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12 .  Thermosonde  17 


Figure  36.  Cl  profile  comparison  for  24OCT01  04Z. 

This  model  run  reasonably  mimics  the  thermosonde 
profile^  as  seen  in  Figure  36.  The  cause  of  the  one 
anomaly  at  11  km  is  uncertain  when  vertical  profiles  are 
compared^  Figures  160  And  178.  No  cloud  ice  was  present 
in  the  model. 
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13.  Thermosonde  18 


Figure  37.  Cl  profile  comparison  for  25OCT01  02Z. 

The  profiles  for  the  case  shown  in  Figure  37  are 
reasonably  similar.  Although  the  levels  are  not  exactly 
matched^  the  trends  are  represented^  Figures  162  And  179. 
Cloud  ice  is  present  in  the  model  at  11  km. 
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14.  Thermosonde  19 


Figure  38.  Cl  profile  comparison  for  25OCT01  04Z. 

At  11  km  the  thermosonde  indicates  significant  drying 
in  Figure  38^  indicative  of  the  stratospheric  boundary^ 
however  the  model  does  not  depict  this  when  Figures  163  and 
180  are  compared. 
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15.  Thermosonde  21  and  Thermosonde  22 


No  comments  are  provided  for  Thermosondes  21  and  22 
(Figures  39  and  40)  because  model  and  thermosonde  profiles 
are  relatively  consistent  with  one  another.  Also^  no 

cloud  ice  was  modeled  in  either  case. 


Thermosoode :vanfft021  txt  Model  c2001 10251 2(01 4  txl 
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Figure  39.  profile  comparison  for  26OCT01  02Z. 


Figure  4 0 . 


g  profile  comparison  for  26OCT01 


04Z  . 
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C.  QUANTITIVE  EVALUATION  OF  Cl  PROFILES  DERIVED  FROM 
THERMOSONDE  AND  COAMPS  DATA 

The  quantitative  evaluation  of  the  model's  performance 
was  determined  by  taking  the  C\  from  both  the  model  and 
thermosonde  data  used  in  Section  B  and  determining  the 
isoplanatic  angle  ( )  for  adaptive  optics^  which  weighs 

mid-  and  upper-troposheric  C\  more  heavily.  The  6^  is  a 
measurement  of  a  plane  over  which  a  source  of  light  would 
undergo  wavefront  distortions.  The  critical  6^  value  that 
separates  severe  from  non-severe  optical  turbulence  is 

^.bJLlrads,  Any  values  that  fall  below  this  is  considered 
severe  optical  turbulence.  Values  that  fall  above  this 
number  are  considered  non-severe.  Therefore^  the  model 

simulations  giving  larger  values  than  what  was  observed  are 
indicating  a  bias  towards  predicting  weaker  turbulence  that 
what  was  measured.  The  equation  used  to  determine  6^  is 
defined  as  (Frieda  1982) 

H  “|“5 

e^  =  2.9\{27rlXf  \h'cl{h)dh  .  (21) 

0 

For  this  case^  the  X  was  set  to  0.5  /Llm  and  the  integration 

was  evaluated  from  the  surface  to  17.5  km.  Furthermore^ 

the  trapezoidal  integration  process^ 

j  fix)  dx  =  ^  (x.^1  +Xi){fM+f)  12,  (22) 

a 

was  utilized  to  determine  the  integration  approximation. 
Figure  41  shows  the  results  from  the  October  observation 
period  and  Table  4  gives  the  values  used. 
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October  Case 


Figure  41 


0^  Scatter  plot  results  for  COAMPS  and 


thermosonde  C"  values  for  October  2001. 


Thermosonde  Thermosonde  Model 


Number 

^0 

2 

0, 

.  1144 

0, 

.2819 

3 

0, 

.0997 

0, 

.2324 

|5 

0, 

.  0778 

0, 

.  1526| 

6 

0, 

.1513 

0, 

.0882 

P 

0, 

.  1111 

0, 

.  0438| 

10 

0, 

.  0418 

0, 

.2292 

111 

0, 

.1203 

0, 

.  1720| 

12 

0, 

.  1074 

0, 

.  1785 

|13 

0, 

.1189 

0, 

.  1784| 

15 

0, 

.1566 

0, 

.1552 

|16 

0, 

.0969 

0, 

.  1032| 

17 

0, 

.  1068 

0, 

.0976 

|18 

0, 

.  1357 

0, 

.  1216| 

19 

0, 

.2066 

0, 

.1875 

|21 

0, 

.2122 

0, 

.29541 

22  0.2005  0.2818 


Table  3.  Isoplanatic  angle  calculations  for  October. 
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Similarly  this  process  was  applied  to  the  March 
observation  period  and  results  shown  in  Figure  42  with 
corresponding  values  in  Table  5. 


March  Casa 


Figure  42 


0^  Scatter  plot  results  for  COAMPS  and 


thermosonde  C"  values  for  March  2002. 


Thermosonde  Thermosonde  Model 


Number 

^0 

0^  (xl.Oe-004) 

6 

0 . 1757 

0.2632 

7 

0.1301 

0.2759 

|8 

0.0912 

0.29551 

9 

0.0780 

0.0733 

111 

0 . 0341 

0 .04491 

12 

0.0490 

0.0316 

|13 

0.0992 

0.20301 

14 

0.1700 

0.1897 

|15 

0.1406 

0.16381 

16 

0.1656 

0.1298 

|17 

0.0823 

0 . 0974| 

18 

0 . 0844 

0.1501 

119 

0.1017 

0.27951 

20 

0.1260 

0.1756 

|21 

0.0699 

0.13871 

22 

0.1396 

0.2463 

^3 

0 . 0581 

0.22451 

24 

0.0976 

0.2185 

;  4. 

Isoplanatic  angle 

calculations  for  March. 
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From  the  plots  shown  in  Figures  41  and  42^  the 
conclusion  can  be  drawn  that  the  model's  overall  tendency 
is  to  over-forecast  using  the  TKE  Method.  This  may  be 

attributed  to  the  smoothing  caused  by  courser  resolutions 
in  the  upper  atmosphere^  as  the  standard  COAMPS  vertical 
resolution  scheme  does. 

To  test  this  hypothesis^  the  model  was  re-run  with  83 
vertical  levels  for  selected  cases  in  October  and  March. 
For  October^  model  runs  corresponding  to  Thermosondes  5-10 
were  analyzed  to  investigate  the  jet  event  that  took  place 
on  October  21^^.  For  March^  the  model  runs  corresponding 
to  Thermosonde  22^23  and  24  were  used  because  a  strong  jet 
was  present.  Qualitative  results  for  October  and  March 

are  given  in  Figures  43  and  44^  respectively. 

Figure  43  indicates  that  the  increased  resolution  did 
improve  two  of  the  four  forecasts  in  the  October  cases^ 
with  one  showing  no  improvement.  These  results  were 

encouraging^  however^  when  applied  to  the  selected  March 
observation  period  (Figure  44)^  the  forecasts  in  all  three 
cases  deteriorated.  This  indicates  that  there  is  promise 
in  improving  turbulence  forecasting  by  increasing  the 
resolution  but  that  resolution  alone  will  not  provide  a 
complete  solution . 
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October  Case 


Figure  43.  Hi-resolut ion  COAMPS  results  for  October 

2001 . 


Thermosonde  Thermosonde  Model  Hi-Res 


Number  0^  0^  0^  (xl.Oe-004) 
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111 
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0 . 1189 

o 

-u 

00 

15 
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0 . 1552 
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0.0969 

0 . 1032| 

17 

0.1068 
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|18 

0.1357 

0 . 1216| 

19 

0.2066 

0.1875 

|21 

0.2122 

0.29541 

22 

0.2005 

0.2818 

Table  5.  Higher  resolution  isoplanatic  angle  results  for 

October  2001 . 
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March  Case 


Figure  44.  Hi-resolut ion  COAMPS  results  for  March  2001. 


Thermosonde  Thermosonde  Model  Hi-Res 
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24 
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Table  6.  Higher 

resolution 

isoplanatic  angle  results  for 

March  2002 . 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  quantitative  results  for  October  seemed  to  show 
that  there  was  improvement  in  forecasting  turbulence  in  two 
of  the  four  cases^  model  runs  corresponding  to  Thermosondes 
7  and  10;  the  model  Thermosonde  10  showing  drastic 
improvement.  One  case^  corresponding  to  Thermosonde  6, 
showed  no  forecast  improvement  but  did  increase  the  amount 
of  modeled  turbulence.  Finally^  the  model  corresponding 
to  Thermosonde  5,  increased  the  optical  turbulence  and 
decreased  the  forecast  quality. 

The  quantitative  results  for  March  all  indicate  that 
increasing  the  resolution  lead  to  increasing  the  turbulence 
modeled^  deviating  farther  from  the  thermosonde  calculated 
optical  turbulence  than  the  low  resolution  values. 

Considering  the  synoptic  patterns  affecting  each  of 
these  cases^  one  factor  that  may  have  attributed  to  the 
different  results  is  the  type  of  jet  over  the  point  of 
interest.  In  the  October  case^  an  offshoot  of  the  STJ  is 
affecting  the  mesoscale^  while  the  PFJ  is  interacting  with 
the  STJ  over  California  in  the  March  case.  The  different 
driving  mechanisms^  height  levels^  and  interaction  can  all 
be  affecting  the  model  structure. 

Although  these  findings  are  inconclusive^  they  do 
indicate  that  increasing  the  vertical  resolution  is  one 
small  part  of  the  larger  puzzle.  When  weighing  the 
increased  computational  cost  caused  by  increasing  the 
vertical  resolution  verses  the  value  added^  this  study 
seems  to  indicate  that  the  cost  may  be  too  high.  More 
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likely^  improving  parameterization  relationships  may  have  a 
greater  impact  at  this  time. 

Other  areas  that  could  be  studied  include;  the  cloud 
ice  parameterization  and  its  affects  on  the  calculation  of 
optical  turbulence^  deviations  in  the  moisture  content 
being  modeled  at  and  above  the  tropopause  as  compared  to 
actual  measurements^  and  the  eddy  diffusivity  to  altitude 
relationship.  Additionally^  a  manuscript  in  preparation  by 
Miller  and  Walters  due  for  publication  in  2004  suggest  that 
the  TKE-Free  Method  could  yield  better  results  in  a  high- 
shear  environment^  therefore  a  detailed  comparative 
analysis  between  the  TKE  and  TKE  Free  Methods  for  the  same 
time  periods  would  be  useful  and  enlightening. 


66 


LIST  OF  REFERENCES 


Beniston^  1998:  From 
Springer^ 328  pp . 

Bougeault^  P.^  C.  D.  Hui^  B. 
1995:  Investigation  of  seeing 
mesoscale  numerical  simulation^ 


Turbulence  to  Climate . 


Fleury^  and  H.  J.  Laurent^ 
by  means  of  an  atmospheric 
Appl.  Opt ,  f  34^  3481-3488. 


Deardorff^  J.  W.^  1980:  Stratocumulus-capped  mixed  layers 

derived  from  a  three-dimensional  models  Bound.  Layer 
Meteor . ,18^  495-527. 


Ellrod^  G.^  cited  2004:  Clear  Air  Turbulence.  [Available 
at  http : //orbit -net . nesdis . noaa . gov/arad/fpdt /CAT . html . ] 


Frieda  D.  L.^  1982:  Anisoplanat ism  in  adaptive  optics^  J. 

Opt.  Soc.  Amer.  A. ^  12,  52-61. 

Gal-Chen^  T.  and  R.  C.  J.  Somerville^  1975:  On  the  use  of  a 
coordinate  transformation  for  the  solution  of  the  Navier- 
Stokes  equations^  J.  Computational  Phys . ,  11,  209-228. 


Hodur^  R.  M.^  1997:  The  Naval  Research  Laboratory's  Coupled 
Ocean/Atmosphere  Mesoscale  Prediction  System  (COAMPS) ^  Mon. 
Wea.  Rev.,  125,  1414-1430. 


67 


Kolmogorov^  A.  N.^  1941:  The  local  structure  of  turbulence 
in  incompressible  viscous  fluid  for  very  large  Reynolds 
numbers^  Doklady  Akad.  Nauk.  USSR^  30^  301. 

Kondo^  0.  Kanechika^  and  N.  Yasuda^  1978:  Heat  and 
momentum  transfers  under  strong  stability  in  the 
atmospheric  surface  layer^  J.  Atmos.  Scl.^  35^  1012-1021. 


Klemp^  J.  B.^  and  R.  B.  Wilhelmson^  1978:  The  Simulation  of 
Three-Dimensional  Convective  Storm  Dynamics^  J.  Atmos. 
Sci.,  35,  1070-1096. 

Kyle,  T.G.,  1993:  Turbulence  and  Optics.  Atmospheric 
Transmission^  Emission  and  Scattering^  Los  Alamos  National 
Laboratory,  Butterworth-Heinemann,  140-155 . 


Lee,  J.  T.,  and  R.  Crane,  1994:  Aeronautical  Meteorology. 
McGraw-Hill  Encyclopedia  of  Science  and  Technology,  8th 
Ed. ,  151-154 . 

Mahalov,  A.,  and  B.  Nicolaenko,  cited  2004:  Clear  air  and 
optical  turbulence  in  a  jet  stream  in  the  airborne  laser 
context . 

[http : //www . nave . hpc . mil/Navigator/Spring03_Feature3 . html . ] 


Mellor,  G.  L.,  1973:  Analytic  prediction  of  stratified 

planetary  surface  layers,  J.  Atmos.  Sci.^  30,  1061-1069. 


68 


_ f  and  T.  Yamada^  1974:  A  hierarchy  of 

turbulence  closure  models  for  planetary  boundary  layers^  J. 
Atmos.  Sci.,  31,  1791-1292. 

_ ,  and  T.  Yamada,  1982:  Development  of  a 

turbulent  closure  model  for  geophysical  fluid  problems, ' 
Rev.  Geophysics  and  Space  Physics^  20,  851-875. 

Meteorological  Office  College,  1996:  Turbulence  and  Gust. 
Source  Book  to  the  Forecasters''  ReferencSf  The  Met.  Office, 
1-12  . 


Miller,  D.  K.,  D.  L.  Walters,  A.  Slaviv,  2001:  Evaluation 
of  turbulent  mixing  length  parameterization  applied  to  the 
case  of  an  approaching  upper-troposheric  trough. 
Preprints  of  the  Ninth  Conference  on  Mesoscale  Processes, 
Ft.  Lauderdale,  FL. 


Nuss,  W.  A.,  and  D.  F.  Titley,  1994:  Use  of  multiquadric 
interpolation  for  meteorological  objective  analysis,  Mon. 
Wea.  Rev.,  122,  1611-1632. 

Rosendal,  H.,  cited  2004:  Clear  air  turbulence  and  the  cut¬ 
off  low .  [http : / /members . aol . com/Rosendalhe/ clearair . htm] 


69 


Tjernstrom^  1993:  Turbulence  length  scales  in  stably 

stratified  free  shear  flow  analyzed  from  slant  aircraft 
profiles.  J.  Appl .  Meteor.^  32,  948-963. 

Walters^  D.  L.^  1995:  Measurement  of  optical  turbulence 

with  higher-order  structure  functions^  Appl.  Opt .  3A, 

1591-1597 . 


_ ,  and  L.W.  Bradford^  1997:  Measurements  of  ro . 

and  6^  :  two  decades  and  18  sites^  Appl.  Opt .  ^  3S,  7876- 

7886. 


_ ,  and  D.  K.  Miller^  1999:  Evolution  of  an 

upper-tropospheric  turbulence  event  -  Comparison  of 
observations  to  numerical  simulations.  Preprints  of  the  13^^ 
Symposium  on  Boundary  Layers  and  Turbulence^  Dallas^  TX. 


Warner^  T.  T.^  Peterson^  R.  A.,  and  Treadon^  R.  E.^  1997:  A 
tutorial  on  lateral  boundary  conditions  as  a  basic  and 
potentially  serious  limitation  to  regional  numerical 
weather  prediction.  B.  of  the  Amer.  Meteor.  Soc.,  18,  2599- 
2617  . 


Wolfe^  W.  F.^  and  G.  J.  Zissis^  1978.  Propagation  Through 
Atmospheric  Turbulence.  The  Infrared  Handbook,  Office  of 
Naval  Research^  Department  of  the  Navy^  Washington  ,  DC^  1- 
49. 


70 


APPENDIX  A.  NCEP  CHARTS 


NCEP/MCP  charts  shown  in  this  section  were  taken  from 
the  Pacific  Region  Synoptic  analysis  archived  by  Naval 
Postgraduate  School's  Meteorology  Department.  It  was 

included  to  provide  OOZ  and  12Z  analysis  for  the  period 
from  18OCT01  to  25OCT01  and  01MAR02  to  07MAR02.  Where 

this  was  not  possible,  the  closest  archived  analysis  was 
used . 

1.  October  2001 


Figure  45. 


NCEP/MPC  OOZ 


Surface  Analysis, 


18OCT01 . 
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Figure  46. 


NCEP/MPC  06Z 


Surface  Analysis, 


18OCT01 . 


Figure  4 7 . 


NCEP/MPC  OOZ 


Surface  Analysis, 


19OCT01 . 
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Figure  48. 


NCEP/MPC  12Z 


Surface  Analysis, 


19OCT01 . 


Figure  49. 


NCEP/MPC  OOZ 


Surface  Analysis, 


20OCT01 . 
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Figure  50. 


NCEP/MPC  12Z 


Surface  Analysis, 


200CT01 . 


Figure  51 . 


NCEP/MPC  OOZ 


Surface  Analysis, 


21OCT01 . 
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Figure  52.  NCEP/MPC  12Z  Surface  Analysis,  21OCT01. 


Figure  53. 


NCEP/MPC  OOZ 


Surface  Analysis, 


22OCT01 . 
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Figure  54.  NCEP/MPC  12Z  Surface  Analysis,  22OCT01. 


Figure  55. 


NCEP/MPC  OOZ 


Surface  Analysis, 


23OCT01 . 
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Figure  56.  NCEP/MPC  06Z  Surface  Analysis,  23OCT01. 


Figure  57 . 


NCEP/MPC  OOZ 


Surface  Analysis, 


24OCT01 . 
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Figure  58. 


NCEP/MPC  06Z 


Surface  Analysis, 


24OCT01 . 


Figure  59. 


NCEP/MPC  OOZ 


Surface  Analysis, 


25OCT01 . 
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Figure  60. 


NCEP/MPC  12Z 


Surface  Analysis, 


25OCT01 . 
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2 


March  2002 


Figure  61. 


NCEP/MPC  06Z 


Surface  Analysis, 


01MAR02 . 


Figure  62 . 


NCEP/MPC  12Z 


Surface  Analysis, 


01MAR02 . 
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Figure  63. 


NCEP/MPC  06Z 


Surface  Analysis, 


02MAR02 . 


Figure  64 . 


NCEP/MPC  12Z 


Surface  Analysis, 


02MAR02 . 


81 


Figure  65. 


NCEP/MPC  06Z 


Surface  Analysis, 


03MAR02 . 


Figure  66. 


NCEP/MPC  OOZ 


Surface  Analysis, 


04MAR02 . 
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Figure  67 . 


NCEP/MPC  12Z 


Surface  Analysis, 


04MAR02 . 


Figure  68. 


NCEP/MPC  OOZ 


Surface  Analysis, 


05MAR02 . 
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Figure  69. 


NCEP/MPC  12Z 


Surface  Analysis, 


05MAR02 . 


Figure  7 0 . 


NCEP/MPC  OOZ 


Surface  Analysis, 


06MAR02 . 
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Figure  71.  NCEP/MPC  12Z  Surface  Analysis,  06MAR02. 


Figure  72 . 


NCEP/MPC  OOZ 


Surface  Analysis, 


07MAR02 . 
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APPENDIX  B.  SATELLITE  IMAGERY 


Figure  73.  GOES  IR  Composite  OOZ,  18OCT01. 


Figure  74.  GOES  IR  Composite  12Z,  18OCT01. 
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Figure  75.  GOES  IR  Composite  OOZ,  19OCT01. 


Figure  76.  GOES  IR  Composite  12Z,  19OCT01. 
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Figure  77.  GOES  IR  Composite  OOZ,  20OCT01. 


Figure  7 8 . 


GOES  IR  Composite  12Z,  20OCT01. 
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Figure  79.  GOES  IR  Composite  OOZ,  21OCT01. 


Figure  80.  GOES  IR  Composite  12Z,  21OCT01. 
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Figure  81.  GOES  IR  Composite  OOZ,  22OCT01. 


Figure  82.  GOES  IR  Composite  12Z,  22OCT01. 
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Figure  83.  GOES  IR  Composite  OOZ,  23OCT01. 


Figure  84.  GOES  IR  Composite  12Z,  23OCT01. 
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Figure  85.  GOES  IR  Composite  OOZ,  24OCT01. 


Figure  86.  GOES  IR  Composite  12Z,  24OCT01. 
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Figure  87.  GOES  IR  Composite  OOZ,  25OCT01. 


Figure  88 . 


GOES  IR  Composite  12Z,  25OCT01. 
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APPENDIX  C.  COAMPS  81  KM  RESOLUTION  MODEL 


95 


Figure  94 . 


21OCT01 . 
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Figure  95.  COAMPS  analysis  valid  for  12Z,  22OCT01. 


Figure  96. 
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valid  for  OOZ,  25OCT01. 
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valid  for  12Z,  25OCT01. 
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APPENDIX  D.  AVN  JET  LEVEL  ANALYSIS 


Figure  102. 


AVN  Jet  Stream  Analysis  OOZ,  18OCT01 


Figure  103. 


AVN  Jet  Stream  Analysis  12Z, 


18OCT01 
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Figure  104.  AVN  Jet  Stream  Analysis  OOZ,  19OCT01 


Figure  105. 


AVN  Jet  Stream  Analysis  12Z, 


19OCT01 
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Figure  106.  AVN  Jet  Stream  Analysis  OOZ,  20OCT01 


Figure  107 . 


AVN  Jet  Stream  Analysis  12Z, 


20OCT01 
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Figure  108.  AVN  Jet  Stream  Analysis  OOZ,  21OCT01 


Figure  109. 


AVN  Jet  Stream  Analysis  12Z, 


21OCT01 
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Figure  110.  AVN  Jet  Stream  Analysis  OOZ,  22OCT01 


Figure  111. 


AVN  Jet 


Stream  Analysis 


12Z, 


22OCT01 
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Figure  112.  AVN  Jet  Stream  Analysis  OOZ,  23OCT0 


Figure  113. 


AVN  Jet 


Stream  Analysis 


12Z, 


23OCT0 
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Figure  114.  AVN  Jet  Stream  Analysis  OOZ,  24OCT01 


Figure  115. 


AVN  Jet  Stream  Analysis  12Z, 


24OCT01 


107 


Figure  116.  AVN  Jet  Stream  Analysis  OOZ,  25OCT01 


Figure  117 . 


AVN  Jet  Stream  Analysis  12Z, 


25OCT01 
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APPENDIX  E.  COAMPS  JET  LEVEL  MODEL 


(right)  valid  for  OOZ,  19OCT01. 


(right)  valid  for  12Z,  19OCT01. 


(right)  valid  for  OOZ,  20OCT01. 
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(right)  valid  for  12Z,  200CT01. 


(right)  valid  for  OOZ,  21OCT01. 
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Figure  124.  COAMPS  300  mb  analysis  valid  12Z  22OCT01. 


(right)  valid  for  OOZ,  23OCT01. 


(right)  valid  for  12Z,  23OCT01. 
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(right)  valid  for  OOZ,  24OCT01. 


(right)  valid  for  12Z,  24OCT01. 


(right)  valid  for  OOZ,  25OCT01. 
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Figure  130.  COAMPS  300  mb  analysis  valid  12Z  25OCT01. 
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APPENDIX  F.  RAWINSONDE  PROFILES 
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Figure  131.  Vandenberg  AFB,  18Z  18OCT01. 
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Vandenberg  AFB,  12Z  19OCT01. 
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Figure  133 
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Vandenberg  AFB,  12Z  21OCT01. 
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Figure  139.  Vandenberg  AFB^  12Z  22OCT01. 
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Vandenberg  AFB,  OOZ  23OCT01. 
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LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  280.2 
LCLP  865.9 
MLTH  292.0 
MLMR  7.41 
THCK  5711. 
PWAT  23.71 


12Z  23  Oct  2001 


University  of  Wyoming 


Figure  141.  Vandenberg  AFB^  12Z  23OCT01. 


SLAT  34.75 
SLON  -120. 
SELV  121.0 
SHOW  8.27 
LIFT  13.17 
LFTV  13.15 
SWET  112.1 
KINX  8.10 
CTOT  12.10 
VTOT  23.10 
TOTL  35.20 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  -9999 
EQTV  -9999 
LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  281.5 
LCLP  905.0 
MLTH  289.7 
MLMR  7.79 
THCK  5722. 
PWAT  21.60 
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Figure  142 . 


Vandenberg  AFB,  OOZ  24OCT01. 
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SLAT  34.75 
SLON  -120. 
SELV  121.0 
SHOW  10.27 
LIFT  1 5.44 
LFTV  15.39 
SWET  35.94 
KINX  -4.10 
CTOT  2.90 
VTOT  27.90 
TOTL  30.80 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  -9999 
EQTV  -9999 
LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  271.9 
LCLP  779.8 
MLTH  292.0 
MLMR  4.54 
THCK  5722. 
PWAT  12.91 
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Figure  143.  Vandenberg  AFB^  12Z  24OCT01. 


72393  VBG  Vandenberg  Afb 


SLAT  34.75 
SLON  -120. 
SELV  121.0 
SHOW  4.57 
LIFT  8.87 
LFTV  8.79 
SWET  33.84 
KINX  13.70 
CTOT  13.90 
VTOT  28.90 
TOTL  42.80 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  -9999 
EQTV  -9999 
LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  279.6 
LCLP  862.5 
MLTH  291.7 
MLMR  7.16 
THCK  5695. 
PWAT  20.87 
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Figure  144 . 


Vandenberg  AFB,  12Z  25OCT01. 
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APPENDIX  G.  COAMPS  VERTICAL  PROFILE  FORECASTS 


Figure  145.  COAMPS  forecast  for  OOZ,  19OCT01. 


Figure  146.  COAMPS  forecast  for  02Z,  19OCT01. 
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Figure  147.  COAMPS  forecast  for  04Z,  19OCT01. 


Figure  148.  COAMPS  forecast  for  OOZ,  20OCT01. 
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Figure  149.  COAMPS  forecast  for  02Z,  200CT01 


Figure  150.  COAMPS  forecast  for  04Z,  200CT01 


125 


Figure  151 . 


COAMPS  forecast  for  05Z,  200CT01 


Figure  152.  COAMPS  forecast  for  OOZ,  21OCT01 
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Figure  153.  COAMPS  forecast  for  02Z,  210C 


Figure  154.  COAMPS  forecast  for  04Z,  21OCT01 
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Figure 


155.  COAMPS  forecast  for  OOZ,  23OCT01. 


Figure  156.  COAMPS  forecast  for  02Z,  23OCT01. 
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Figure  157.  COAMPS  forecast  for  04Z,  23OCT01 


Figure  158.  COAMPS  forecast  for  OOZ,  24OCT01 
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Figure  159.  COAMPS  forecast  for  02Z,  24OCT01. 


Figure  160.  COAMPS  forecast  for  04Z,  24OCT01. 
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Figure  161.  COAMPS  forecast  for  OOZ,  25OCT01. 


Figure  162.  COAMPS  forecast  for  02Z,  25OCT01. 
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Figure  163.  COAMPS  forecast  for  04Z,  25OCT01. 


Figure  164.  COAMPS  forecast  for  OOZ,  26OCT01. 
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Figure  165.  COAMPS  forecast  for  02Z,  26OCT01. 


Figure  166.  COAMPS  forecast  for  04Z,  26OCT01. 
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APPENDIX  H.  THERMOSONDE  PROFILES 


Figure 


20041002/00  74606  VBG  SELV:  112  LIFT:  11 KINX:  12TOTL:  35 
SWET:  28LCLP:  976  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


167.  Thermosonde  2,  launch  time  0132Z  19OCT01 


20041003/00  74606  VBG  SELV:  112  LIFT:  11  KINX:  10  TOIL:  33 
SWET:  35LCLP:  976  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 
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20041005/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  10  TOIL:  36 
SWET:  37LCLP:  973  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  169.  Thermosonde  5^  launch  time  0115Z  200CT01 


20041006/00  74606  VBG  SELV:  112  LIFT:  8  KINX:  13TOTL:  39 
SWET:  31LCLP:  973  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  170.  Thermosonde  6^  launch  time  0300Z  200CT01 
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20041007/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  9TOTL:  37 
SWET:  -9999  LCLP:  974  LFCV:  -9999  EQLV:  -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  171.  Thermosonde  1,  launch  time  0444Z  200CT01. 


20041010/00  74606  VBG  SELV:  112  LIFT:  12  KINX:  21TOTL: 

SWET:  82  LCLP:  967  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


42 


Figure  172.  Thermosonde  10^  launch  time  0235Z  21OCT01. 
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Figure 


20041011/00  74606  VBG  SELV:  112  LIFT:  13  KINX:  24  TOIL: 

SWET:  107LCLP:  968  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


42 


173.  Thermosonde  11^  launch  time  0413Z  21OCT01 


Figure 


20041012/00  74606  VBG  SELV:  112  LIFT:  12  KINX:  -9TOTL: 

SWET:  -9999  LCLP:  967  LFCV:  -9999  EQLV:  -9999 
CAPE:  0  CINS:  0  BRCH:  0 
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174.  Thermosonde  12^  launch  time  0115Z  23OCT01 
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20041013/00  74606  VBG  SELV:  112  LIFT:  11  KINX:  10  TOIL: 

SWET:  64LCLP:  970  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


32 


Figure 


20041015/00  74606  VBG  SELV:  112  LIFT:  13  KINX:  -18TOTL:  23 
SWET:  69LCLP:  948  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


176.  Thermosonde  15^  launch  time  0114Z  24OCT01 
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20041016/00  74606  VBG  SELV:  112  LIFT:  12  KINX:  -22TOTL:  15 
SWET:  75LCLP:  954  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  177.  Thermosonde  16^  launch  time  0248Z  24OCT01 


20041017/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  -22TOTL:  17 
SWET:  61LCLP:  966  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  178.  Thermosonde  11,  launch  time  0430Z  24OCT01 
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20041018/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  -7TOTL:  37 
SWET:  -9999  LCLP:  873  LFCV:  -9999  EQLV:  -9999 
CAPE:  0  CINS:  0  BRCH:  0 
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Figure  179. 


Thermosonde  18^  launch  time  0115Z  25OCT01 


20041019/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  -IITOTL:  33 
SWET:  40  LCLP:  903  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 
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Figure  180 . 


Thermosonde  19,  launch  time  0253Z  25OCT01 
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20041021/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  7TOTL:  34 
SWET:  33LCLP:  979  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  181.  Thermosonde  21,  launch  time  0249Z  26OCT01 


20041022/00  74606  VBG  SELV:  112  LIFT:  10  KINX:  6TOTL:  33 
SWET:  34LCLP:  977  LFCV: -9999  EQLV: -9999 
CAPE:  0  CINS:  0  BRCH:  0 


Figure  182.  Thermosonde  22^  launch  time  0423Z  26OCT01 
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